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The three-dimensional nature of plasmoid instabilities is studied using the reduced magnetohydrodynamic 
equations. For a Harris equilibrium with guide field, represented by = Bpotdinh{x/X)y-\- BzoZ^ a spectrum 
of modes are unstable at multiple resonant surfaces in the current sheet, rather than just the null surface of 
the polodial field Byo{x) = tanh(x/A), which is the only resonant surface in 2D or in the absence of a 
guide field. Here Bpo is the asymptotic value of the equilibrium poloidal field, Bzo is the constant equilibrium 
guide field, and A is the current sheet width. Plasmoids on each resonant surface have a unique angle of 
obliquity = diTctdin{kz/ky). The resonant surface location for angle 6 is Xs = —Aarctanh(tan 6>52o/^po), and 
the existence of a resonant surface requires |^| < diTctdin{Bpo/ B^o)- The most unstable angle is oblique, i.e. 
0^0 and Xg ^ 0, in the constant-?/^ regime, but parallel, i.e. ^ = and Xs = 0, in the nonconstant-?/^ regime. 
For a fixed angle of obliquity, the most unstable wavenumber lies at the intersection of the constant-?/; and 
nonconstant-?/; regimes. The growth rate of this mode is 7max/ro — ^^^"^{1 — /i^)^/^, in which To = Va/L^ Va 
is the Alfven speed, L is the current sheet length, and Sl is the Lundquist number. The number of plasmoids 
scales as - 5^/^(1 - + /i')^/^ 

PACS numbers: 52.35.Vd,52.22.Tn,94.30.cp,96.60.Iv 



I. INTRODUCTION 

Plasmoid dominated reconnection occurs when a thin 
current sheet breaks into a chain of secondary islands, 
or plasmoids, which convect along the reconnection out- 
fiow, eventually removing magnetic field from the cur- 
rent sheet .^^^^ Before being ejected, plasmoids may 
coalesce,!^ and current sheets between plasmoids can 
excite new generations of plasmoids. 1^^^ This process 
may be con sidere d turbulent if it can repeat sufficiently 
many times .f^^ESI Since plasmoid dominated reconnec- 
tion proceeds much faster than the conventional Sweet- 
Parkeill^^ rate, onset of the instability can trigger fast 
reconnection. In 3D plasmoids are tube-like in shape, see 
Fig.[l] and are often called fiux ropes. There is significant 
observational evidence for fiux ropes in astrophysical re- 
connection sites including solar fiares,1^ and the Earth's 
magnetopause^^ and magnet ot ail. Similar tearing in- 
stabilities are also important in magnetic confinement 
fusion experiments. 

Magnetohydrodynamic (MHD) theories of tearing in- 
stabilities are typically 2D in that they take the guide 
field direction (i) to be ignorable.-^^^ For these parallel 
modes, kz = 0, fiux ropes are aligned with the guide field 
as shown in Fig.[l^. Oblique tearing modes, with kz ^ 0, 
are a 3D effect. In tokamak parlance, oblique modes are 
those with n 7^ 0, where n is the toroidal mode num- 
ber. The primary differences between tokamak tearing 
modes'^ and plasmoid instabilities^ ^ are the current dis- 
tributions and boundary conditions. Tearing modes are 
instabilities of diffuse current distributions, whereas plas- 
moids are secondary instabilities of thin current sheets. 
The difference has consequences for how the instabilities 
scale with resistivity. Boundaries are periodic in both the 
toroidal and poloidal directions in a tokamak, whereas 




FIG. 1. Schematic depiction of constant flux surfaces of plas- 
moids in 3D. Dashed lines designate the guide field direction 
(z). Flux ropes for parallel modes (a) are aligned with the 
guide field, while for oblique modes (b) they are misaligned 
by angle 0. 



in astrophysical situations in which current sheets arise 
boundaries are often open or line-tied. 

It has recently been shown^ that the linear plasmoid 
instabilit}'^ can be related to the conventional tearing 
mode^^ in a simple way. For a Harris equilibrium^^ with- 
out guide field, = 5^tanh(x/A)^, the classical tear- 
ing mode growth rate is^ 



for the constant-?/^ and nonconstant-?/^ regimes, respec- 
tively. The maximum growth rate occurs at the in- 
tersection of the two branches, kXS^^^ o:^ 1, where 
Tmax^A ~ 6'""^/^. Here, S = tr/tais the Lundquist num- 
ber based on the current sheet width, tr = AttX'^ /{rjc^) 
is the resistive diffusion time, ta = X/Va = XyJ A-k p / Bpo 
is the Alfven time, and k is the wavenumber. The re- 
cent insight connecting the plasmoid growth rat^ with 
Eq. was to account for the Lundquist number seal- 
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ing of the current sheet widthP MHD current sheets 

— 1/2 

obey the Sweet-Parker width A = ^sp = LS^^ ' , in 
which Sl = ^"kLVa/ {c^Vi) = {L/X)S is the Lundquist 
number based on the current sheet length. For a cur- 
rent sheet, the maximum growth rate from Eq. ([T]) is 

7max — S\/^Va/L^ which scales with S'^ to a positive 
exponent, rather than the negative exponent scaling of 
the most unstable tearing mode. Because current sheets 
becomes increasingly singular at high Sl^ the plasmoid 
growth rate is large for high Sl plasmas such as the solar 
corona {Sl ^ 10^^) and fusion experiments {Sl ^ 10^). 
Plasmoid growth rates that scale with Sl to a positive 
exponent have also been calculated in the Hall MHD 
regime. 

The present work is motivated by a recent study by 
Daughton et al.,^^ who showed using linear Vlasov and 
particle-in-cell simulations of the Harris current sheet^^ 
that the conventional collisionless kinetic theorjEl breaks 
down for oblique modes. They suggest that this is a 
failure of asymptotic boundary layer analysis, which is 
a challenge to analytic theory. Here, we consider the 
simpler resistive MHD problem and show that, at least 
within this framework, boundary layer theory can faith- 
fully describe oblique tearing modes. This regime is im- 
portant in its own right since many reconnection prob- 
lems of interest are sufficiently collisional that a resis- 
tive MHD model is warranted. The spectrum of oblique 
tearing modes has important consequences for the gen- 
eration of turbulence by overlapping flux ropes It is 
also important when considering whether plasmoids can 
fill the volume of a current sheet, which is an important 
assumption in some particle acceleration theories.^ We 
find that plasmoids are volume filling, but the angle of 
obliquity and growth rate depend on the resonant surface 
location. A numerical study of oblique tearing using Hall 
MHD without guide field has been presented by Cao and 
Kan.'^ Huang and Zweibel^^ studies the reduced MHD 
problem numerically, with guide field, focusing on line- 
tied boundary conditions. 

Plasmoid instabilities arise at resonant surfaces of the 
ideal MHD equations, defined by k • = 0, near 
which dissipation (resistivity here) allows for reconnec- 
tion of magnetic field. ^ The primary difference be- 
tween the 2D and 3D theories is the location of res- 
onant surfaces. For a Harris sheet with guide field. 
Bo = Bpo tdinh{x / X)y + BzoZ^ resonant surfaces are lo- 
cated at Xs = — Aarctanh(tan^52o/5po), where 6 = 
diTctdiii{kz/ky) is the angle of obliquity. In the conven- 
tional 2D theories, either Bzo = or /c^ = 0, in which 
case there is a single resonant surface corresponding to 
the null of the sheared field Xs = 0. In the 3D problem 
resonant surfaces can be found at any location across 
the current sheet. Modes at each surface correspond to 
different angles of obliquity. The angle for a mode at 
surface location x is 6 = arctan[tanh(— x/Aj^po/^^o]. 
Parallel modes {6 = 0) are found at the null surface 
of the poloidal field and the magnitude of the angle 
of the mode increases with distance from the null sur- 



face. The existence of a resonant surface requires \0\ < 
diTctdin{Bpo/Bzo)' We show that for large k (the constant- 
tjj regime) the most unstable mode is oblique, satisfying 
:L{Bpo/ Bzo)\/{'^ + /c^A^)/3, and that parallel modes 
are a local minimum of the growth rate. For small k (the 
nonconstant-?/^ regime), parallel modes are the most un- 
stable, and the growth rate falls off monotonically with 

1^1- 

The rest of this paper is organized as follows. Sec- 
tion IH] describes the reduced MHD equations that form 



the basis of this analysis. Section |III| provides a bound- 
ary layer theory for the tearing mode growth rate, and 
this is used to derive the dispersion relation for the plas- 
moid instability. These results are compared with direct 
numerical solutions of the linearized reduced MHD equa- 
tions in Sec. HVl Section [Vl discusses numerical solutions 
of the flux and stream functions, which change with the 
angle of the mode. Section [Vl| provides a summary of the 
results. 



II. REDUCED MHD EQUATIONS 

The reduced MHD equations are based on tokamak 
ordering^^ 



d:,,dy,Bzo-0{l), 

dz,dt,^l)o,(l)o,Vzo ^ 0{e), 
BzuVzi,^i,^i^O{e^), 



(2) 



assuming a constant plasma density, a strong, constant, 
guide field (in the z direction here), and that wavelengths 
in the guide field direction are much longer than in the 
perpendicular directions. 

Applying these approximations, the MHD equation of 
motion 



Ohm's law 



(S, + V • V) V = J X B - VP, 



E + V X B = S-^J, 



(3) 



(4) 



and the relevant Maxwell equations V-B = 0, VxE = 
-dtB, and V X B = J, lead to the reduced MHD 
equations: 
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dtn^in,^] = [Jz,^]^BzdzJz. (5) 



dt^ = Bzdz^ + [(/>, ^] + S-^\/l^ + Eo. (6) 

Here, the stream function is defined by V = x 
z + VzZ^ the flux function by B = V^V^ x i + BzZ, 
Q = — is the vorticity, Jz = — V^V^ is the elec- 
tric current in the z direction, Eo and Vq are constants 
of integration, = dxX + dyy is the perpendicular gra- 
dient, and [f^g] = (V/ x V^) • z is the Poisson bracket. 
Spatial scales are normalized to the current sheet width 
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(x = x/A), velocities to the Alfven speed (V = V/Va), 
time to the Alfven time i = Va^I\ magnetic field to 
the magnitude of the asymptotic poloidal magnetic field 
(B = B/5po), and currents by J = J/[c5po/(47rA)]. 
Tildes have been omitted for notational convenience. 

We linearize Eqs. ([5| and ^ according to f = fo-\-Sf 
in which ipo = i^oi^), (j^o = <l^o(x^ y), and B^o is constant. 
We also assume that fiow profiles satisfy, V^(/)o = 0, and 
that the instability growth rate is much larger than the 
timescale for equilibrium fiows, dxcj^o^ dy(j)o <C 7. Pertur- 
bations satisfy Sf = fi{x) ex.p[i{kyy + kzz) + jt]. Ap- 
plying this procedure, and the tokamak ordering Q, the 
linearized form of Eqs. ([5| and ([6| are: 

j{4>'l -kl4>i) = iF{4>'{ -kl4>,)- iF"^i , (7) 




-iiPi=iF4>i+S-\^'(-kl^^), 



(8) 



in which F = k • and primes denote x derivatives. 
Tokamak ordering implies k = \ky\[l-\-0{e)] and kz/ky = 
tan(6>) ~ 6> + O(e^), where ~ 0{e). In the remainder 
of this work, we adopt the notation {k,6) in place of 
{ky, kz). Thus, F = k Bo = k[Boy{x) + OBoz], which is 
an 0{e) quantity. 



III. BOUNDARY LAYER ANALYSIS 

A. Outer Region 

In the outer region, we assume <C 7 <C 1, in which 
case Eq. ^ reduces to the classical ideal MHD outer 
region of Furth, Killeen, and Rosenbluth (FKR)^ 



(9) 



Equation ([9| holds everywhere except a small region 
about the resonant surface, where F = 0. We follow 
the conventional boundary layer analysis, which matches 
the jump in the first derivative of denoted the tearing 
stability index 



A'^[i>[{xt)-i>[{x;)]/Mxs), 



(10) 



in the inner and outer regions. Here xf = lim^^o{xs^e) 
and ipi is continuous at Xg. 

FKR provides an asymptotic analysis for A' in the 



For k^ > di 



so 



ipi{xs) exp(— 

For P <^ dl, FKR show 



the solu 



large and small k limits, 
tion of Eq. ([9| is = 
A' -2k in this limit. 
A' {l/k)[F'{xs)]'^{FZ^^F-^). An approximate solu- 
tion that captures both the large and small k limits can 
be obtained by adding the asymptotic solutions 



A'^(aVfc)(F_-^ + F-2)-2fc. (11) 

in which a = F' {x = Xs) = kB'^y^x — Xg). 

For the Harris equilibrium with a guide field, B^ = 



5potanh(x) + Bzo^ and Xs 



-Aarctanh(/i) where /i 



FIG. 2. (Color online) Angular dependence of the tearing 
stability index for a Harris equilibrium with Bpo/Bzo = 0.1, 
and four different wavelengths: k\ — 0.4,0.7,1.2 and 2.0. 
Black circles show a numerical solution of Eq. (|9]), blue dashed 
lines the FKR approximation from Eq. ( |12| ), and red dash- 
dotted lines Daughton's approximation from Eq. (|13|). 
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FIG. 3. (Color online) Wavenumber dependance of the tear- 
ing stability index for a Harris equilibrium with Bpo/Bzo = 
0.1, and three values of the angle of obliquity: — 0.00,0.05 
and 0.09 radians. Black circles show a numerical solution 
of Eq. (9l, blue dashed lines the FKR approximation from 
Eq. ([12), and red dash-dotted lines Daughton's approxima- 
tion from Eq. ( 13 ). 



k,B,o/{kAo) - OB,o/Bpo. Thus, a = kBp,{l - /x^) 



and Eq. yields 



AiiC^2[(l + M')/fc-fc]. 



(12) 



Figures [2] and [s] show that Eq. ([T2| agrees well with 
numerical solutions of Eq. ([9| for the Harris equilib- 
rium. Solutions are shown for Bpo/Bzo = 0.1, at fixed 
wavenumbers varying the angle (Fig.[2|, and at fixed an- 
gles varying the wavenumber (Fig.p|. Daughton et al^ 
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have also proposed the solution 



A' 



:( l-fc/2) 
1-k 



(13) 



for a Harris equilibrium with guide field. Predictions of 

Handjs] 

Eqs. (|12|) and (13|) both asymptote to 



Eq. (13) are also shown in Figs. 
For small /c, 

A' 2(1 + However, foF" large k, ^ -2A:, 

while AJ^ ^ — /c(2 + /i^). The /i dependence of the large 
k limit of Eq. lis]) is incorrect, as the asymptotic so- 
lution of Eq. Q and Fig. [s] show. However, Eq. (13) 
provides an excellent approximation for small k. This is 
typically the most interesting case since tearing instabil- 
ity requires A' > and A' becomes negative for large k. 
The simple expression (12) provides an adequate approx- 



imation for all k and ^, capturing both asymptotic limits. 
Both results are exact for normal modes {0 = 0).^^ Fig- 
ure |2] shows that an interesting situation can arise for 
/c ~ 1, where oblique modes are unstable, A\6 ^ 0) > 0, 
but normal modes are stable, A\6 0) < 0. This fea- 
ture is discussed in detail in Sec. |IV| Figures |2]and|3] 
include only \6\ < 0.1, since there is no resonant surface 
for 1^1 >Bpo/B,o = 0.1. 



B. Inner Region 

In the inner region, x — = ^ <C 1, we assume 9^ = 
9| ky and expand F to linear order about the resonant 
surface: F F'{xs){x — Xg) = o^^. Here, Eqs. ^ and 
(|8| reduce to 



and 



7t/-i-ae(i0i) = 5-Vi- 



(14) 



(15) 



Equations ( |14[) and (15) are the same equations used by 
Coppi et a/.^^to calculate A' in the inner layer. We pro- 
vide an alternative derivation using a Fourier transform 
method similar t o tha t developed by Bondeson et al.^^ 
and Porcelli et al^^^ 

The fourth-order system of equations ( 14) and ([l5| has 
solutions with ipi constant, linear, and a solution where 
V^i is an even function of ^ — We are 

interested in the last of these. For large ^, Eqs. ([l4| and 
( 15 ) reduce to 



(16) 



which has the solution = ai^ + arctan(a^/7). For 
large ^, the asymptotic limit of this is 



(17) 



where A and B are constants. Equations (10) and (17) 
imply 



A' = 2A/B. 



The coefficients A and B can be calculated by match- 
ing Eq. (17) with an exact solution of Eqs. (14) and 
(15). It is convenient to do this matching in a Fourier- 



transformed space: f{p) = exp(— Ap- 

plying this, Eq. (17) can be written 



^ + 5(2^J(p) + g|p|)+(9(p^). (19) 



P 



The fourth order system of equations (14) and (15) can 
be written as a second order equation for lipi 



d 

dp 



1 d 
dp 



(20) 



in which A = jS^^^a and p 
solution of Eq. ( 20 ) is 



(S'a)-i/3p. The 



P 



C2 

"^-.3/4 



H",.../4(\/Ap2), 



(21) 



in which M and W are the Whittaker-M and W func- 
tions, u = — A^/^/4, and Ci and C2 are constants. 

The first term of Eq. (21) diverges for large so we 



must take Ci = 0. For p <C 1, 

W^..3/4/p'/' = «(p-' + A^/V2) + b\p\ + 0{p^) (22) 

in which a = ^/^/{2K^'^T[{K^/'^ + 5)/4]} and b = 
40FA5/8/{3r[(A3/2 - l)/4]}. Matching the coefficients 
of Eqs. (|l9j) and (El]) gives A = -C2{Sa)^a/2, and 
B = 3C2(Say/H/(TTA'^). With these, Eq. (isj) provides 
the dispersion relation 



A' 



.l(SaY/SA^/^m^!l^im .23) 



Equation (23) was first obtained by Coppi et al.^^ 



Analytically tractable solutions for the growth rate can 
be obtained from Eq. (23) in the limits A <C 1 (the 
constant-?/^ regime) 



r(i/4) 



4/5 



^-3/5^2/5^/4/5 



and A 



_27rr(3/4)_ 
(the nonconstant-V^ regime) 



3A' 



(24) 



(25) 



The second term in Eq. ([25| causes stabilization at large 
/c, but is typically negligible for the most unstable mode. 



IV. LINEAR GROWTH RATE 

The dispersion relation for the linear tearing mode 
growth rate is obtained by equating Eqs. (12) and (23). 



(18) These equations can also be used to derive the plasmoid 
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dispersion relation for a Sweet-Parker current sheet, as 
was done for parallel modes in Ref. 2^ using the method 
summarized in Sec. [T[ In terms of unnormalized units, 
the tearing mode growth rate in the constant-?/^ and 
nonconstant-?/; regimes from Eqs. (I24|) and (25) are 



5-3/^(i^A)-2/5(l - ^2)2/5(1 ^ ^2 _ ^2 ^2)4/5 



5-l/3(i^A)2/3(l- ^2)2/3^ 



(26) 



For a fixed angle of obliquity, the most unstable 
wavenumber occurs at the intersection of the two 
regimes, which is k\ — + 

^2^3/4 Constant-?/^ corresponds to kX larger than 
this value, and nonconstant-?/; to kX smaller than this 
value. The leading constant coefficients in Eq. (26) are 
[r(l/4)/(7rr(3/4))]^/^ 0.95 for the constant-?/^ regime, 
and unity for the nonconstant-?/; regime. Equation (26) 



reduces to the classical tearing mode dispersion relation 
of Eq. ([T]) for normal modes (/i = 0). The maximum 
growth rate is TmaxTA ^ - ja^y^^. 

Oblique modes (/i / 0) are the most unstable in the 
constant-?/^ regime, where the maximum growth rate has 
angle tan (9 = ±{Bpo/B^o)Vi^ + k'^X^)/3. Here, parallel 
modes are a local minimum in the growth rate. This be- 
havior should be contrasted with that in the nonconstant- 
tjj regime, where parallel modes (/i = 0) are the most 
unstable, and the growth rate falls off monotonically for 
oblique angles. In both regimes, /i < 1 is required for 
instability, otherwise there is no resonant surface. 

The plasmoid growth rate can be calculated from the 
tearing mode dispersion relation by taking the equilib- 
rium to be a Sweet-Parker current sheet, which has width 



A 



ysp 



LS 



-1/2 



L 

1/2. 



With this, 5' = {X/L)Sl = S\ 



1/2 



and ta = 1/(5'^ To), where Tq = Va/L. The plasmoid 
growth rate in the constant-^^ and nonconstant-?/^ regimes 
can then be written 



7/ro 



;2/5^_2/5 
/t2/3(l 



Sr^n-'^/^l - ^2)2/5(1 + ^2 _ ^2/5^)4/5^^^ 



in which n = kL. 

The most unstable angle in the constant-?/; regime is 
e ^ ±{Bpo/B^o)Vi^ + i^^/Sl)/S in these variables. The 
two regimes meet at 



(28) 



where the maximum growth rate is 



7n 



4x1/2 



Equation (28) provides an estimate for the number of 
plasmoids expected to initially arise in an unstable cur- 
rent sheet: ^ /^max/(27r). Parallel modes generate the 
fewest number of plasmoids, and the plasmoid number 
increases monotonically with 0. Equations (28) and (29) 



reduce to the results of Loureiro et al} for parallel modes. 
Equation (27) also provides the instability criterion 



2x1/2 



(30) 
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FIG. 4. (Color online) Wavenumber dependence of the plas- 
moid growth rate for a Harris equilibrium and Bpo/Bzo = 0.1. 
(a) Normal modes {0 — 0) for three values of the Lundquist 
number Sl = 10^, 10^, and 10^°. Circles show the growth 
rate from a direct numerical solution of Eqs. (I7| and ([8|, red 
dashed lines from the boundary layer theory of Eqs. \V2\ and 
( 23 ) , and the blue solid line from boundary layer theory using 
the small k limit of Eq. ^ [A'A ^ 2/(A;A)] and Eq. ([23]). (b) 
Oblique modes with angles = 0.00,0.06,0.09 and 0.099 at 
fixed Sl — 10^. Lines represent the boundary layer theory 
and circles the numerical solutions. 



for the wavenumber. 

Figures |4] and [5] show the plasmoid growth rate calcu- 
lated from the boundary layer theory of Eqs. (12) and 
(23), as well as from a direct numerical solution of the 



(^29) linear reduced MHD equations, ^ and ([8|. All figures 
use Bpo/Bzo = 0.1. Figure shows excellent agree- 



ment between the theory and numerical results for par- 
allel modes at three values of Sl- The k,'^^^ scaling of the 
nonconstant-?/; regime, and k~'^I^ scaling of the small k 
limit {n^ / Sl <C 1) of the constant-?/^ regime are also con- 
firmed. The solid lines in Fig. [4^ are obtained using the 
smah kX limit of Eq. ([l2| [A'A = 2/(A:A)], which extends 
the K~^l^ scaling beyond its region of validity. This is 
the limit assumed in Ref. 1, and is shown for comparison. 
Although the growth rate falls off rapidly for k > /^max, 
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FIG. 5. (Color online) Angular dependence of the plasmoid 
growth rate for a Harris equilibrium and Bpo/Bzo = 0.1. 
Lines represent solutions from the boundary layer theory of 
Eqs. |T2| and j23}, and circles from direct numerical solu- 
tions of Eqs. I?! and (|8|. (a) For fixed Lundquist num- 
ber Sl — 10^, and four values of the wavenumber: n = 
1 X 10^, 5 X 10^, 9 X 10^, and 1 x 10^. (b) For fixed wavenum- 
ber K, — 1 X 10^ and four values of the Lundquist number 
Sl = Ix 10^ 2 X 10^ 1 X 10^ and 1 x 10^°. 



the scaling holds near the intersection with the 

nonconstant-?/^ regime. The maximum growth rate ob- 
tained from the intersection of the nonconstant-^/^ regime 
and the small k limit of the constant-?/; regime provides a 
good approximation of the full analytic, and numerical, 
results. 

Figure [ZJd shows results for fixed 5'^;, = 10^, and dif- 
ferent values of the angle of obliquity. Here, too, the 
growth rate calculated with boundary layer theory com- 
pares well with the numerical results. The agreement 
becomes less favorable for angles very close to the cut- 
off angle |^max| — Bpo/B^o = 0.1. Figure [IJd shows that 
parallel modes are the most unstable in the nonconstant- 
ip regime. The constant-?/^ regime is found n 10^, 
where oblique modes have larger growth rates than par- 
allel modes. The maximum wavenumber for instability 



FIG. 6. (Color online) Contours of the plasmoid growth 
rate calculated from Eqs. (12) and ( [23| as a function of 
wavenumber and angle of obliquity. Here, Sl — 10^ and 
Bpo/Bzo = 0.1. The dashed line shows the estimated bound- 
ary between the constant-'^ and nonconstant-'0 regimes from 
Eq. (I28l. 



from Eq. ( |3Q[ ) shows excellent agreement for all Lundquist 
numbers and angles shown in Fig|4j 

Figure |5] again shows close agreement between bound- 
ary layer theory and numerically calculated growth rates. 
Figure [5^ shows the angular dependence of the growth 
rate for fixed Lundquist number, Sl = 10^, and four 
values of the wavenumber, whereas Fig. ^jp fixes the 
wavenumber, and varies the Lundquist number. In both 
cases, a transition between the two regimes of the in- 
stability are evident. At small a^, or large Sl^ parallel 
modes are most unstable. Here, the growth rate de- 
creases monotonically with \0\ until the stability thresh- 
old at \6\ Bpo/Bzo = 0.1 is reached, which is indica- 
tive of the nonconstant-?/^ regime. Modes at larger a^, or 
smaller Sl^ are most unstable at oblique angles. Here, 
the most unstable angle agrees with th e prediction of the 
constant-?/; regime: :L{Bpo/ Bzo)\/ + / Sl)/'^- 

A contour plot of the growth rate is shown in Fig. |6] for 
a range of wavenumbers near the peak growth rate and 
the entire domain of unstable angles. Here the Lundquist 
number is fixed at Sl = 10^ and the growth rate was cal- 
culated from the boundary layer theory using Eqs. (12) 
and (23). Again, the angular dependence of modes in 



the constant-?/; regime is evident. Here, a linear scale 
has been used for which allows better resolution of the 
constant-?/; regime. The dashed line shows the estimated 
K at the maximum growth rate from Eq. ( [28] ) , which cor- 
responds to the boundary between the constant-?/; (above 
the dashed line) and nonconstant-?/; (below the dashed 
line) regimes. 



Oblique modes in the constant-?/; regime of Eq. (26) 
are analogous to n > 1 tearing modes in a tokamak. In 
fact, all tearing modes in tokamaks are oblique since the 
presence of n = modes would require an infinite safety 
factor (unless m = 0). Furth, Rutherford, and Selberg 
considered n = 1, constant-?/; tearing modes in a periodic 
cylinder. In cylindrical geometry, the linearization be- 
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comes fi{x)ex.p[i{kyy + kz 



/i(r) ex.p[i{kz + mi9)]^ 



0.35 



where 'd is the poloidal direction and k = n/R is quan- 
tized according to the tokamak major radius R. In this 
case, the angle of obhquity is tan^ = kz/ky Tsu/ {Rm)^ 
where is the minor radial location of the resonant sur- 
face. In terms of the safety factor, q = \rBz/{RB^)\, 
q{rs) = m/n^ and 



tan^ 



s/R 



(31) 



The angles of obliquity are small for large aspect ratio 
tokamaks, and increase with the ratio n/m simply due 
to geometry; flux tubes with higher n numbers, or lower 
m numbers, must travel farther poloidally in a toroidal 
transit. The existence of a resonant surface, and the re- 
sulting growth rate, both depend on the current profile, 
which is significantly different in tokamaks than the Har- 
ris equilibrium we have assumed above. The effect of pe- 
riodic boundary conditions on the above analysis would 
simply require quantization of ky and kz according to the 
length of the current sheet L and the domain size in the 
guide field direction. In this case, hv 0^ 2iTm and Fig. |4] 
shows that there are typically thousands of islands in a 
chain for the current sheet lengths of interest. 




V. NUMERICAL EIGENMODE SOLUTIONS 
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Figure [7| shows numerical solutions of the perturbed 
flux and stream functions for fixed 5^ = 1 x 10^, n = Ix 
10^, and four angles of obliquity: = 0,0.025,0.05, and 
0.09. As the angle of obliquity increases, two significant 
changes occur: the eigenfunction center shifts, and sym- 
metries of the eigenfunctions are lost. The shift in eigen- 
function center corresponds to the shift of the resonant 
surface: Xg/X = — arctanh(/i) —diTctdinh{OBzo/Bpo). 
For the plotted angles, 6 = 0.00,0.025,0.05, and 0.09, 
the predicted resonant surface locations are Xs/X = 
0,-0.025,-0.055, and —1.47, respectively. The gradient 
of the perturbed stream function in Fig. [7^ decreases for 
X > Xs^ but grows for x < Xg, as the angle of obliquity 
increases. Likewise, the gradient of the perturbed flux 
function in Fig. ^ decreases for x > Xg^ but remains 
nearly constant on the x < Xs side. 

Figure [S] shows constant flux surfaces, which generate 
flux ropes in 3D, for a parallel mode (6> = in the top 
row) and an oblique mode {6 = 0.06 radians in the bot- 
tom row). Here, the parameters hc = 10^ and Sl = 10^ 
have been chosen. Also shown in the right column are 
2D cuts of the same data, at z = 0, showing magnetic 
islands. The total flux function is used, which is the sum 
of the equilibrium component, and the perturbed compo- 
nent multiplied by a constant amplitude: ip = ipo -\- ai^i- 
For the values in Fig. Isj is normalized to B^qX and a 
is chosen to be 0.01. The equilibrium component, ?/^o, is 
obtained from the definition = V^V^o x ^ + ^20- This 
implies d^Jjo/dx = — tanh(x/A) for the Harris sheet. 



FIG. 7. (Color online) Perturbed stream function 01 (a) and 
flux function ipi (b) for Sl = I x 10^, k = 1 x 10^, and four 
values of the angle of obliquity: = 0,0.025,0.05, and 0.09 
radians. These eigenfunctions were obtained from a direct 
numerical solution of Eqs. ^ and ([8|. 



so 

V^o = -BpoX ln[cosh(x/A)] , (32) 

is used for i/jq. 

The flux ropes shown in Fig. [Sl correspond with the 
qualitative expectation from Fig.^ For parallel modes, 
they are straight and uniform in the z direction, which is 
consistent with taking z to be an ignorable direction in 
the conventional theory. Here, the 2D cuts of magnetic 
islands are symmetric in x about the resonant surface 
X = Xs = 0. For oblique modes, the flux tubes are shifted 
by angle from the normal in the z direction. Here, 
the 2D cuts of magnetic islands do not possess the x 
symmetry of parallel modes, having a slightly shallower 
gradient for x > Xs than for x < Xs- The resonant surface 
is also shifted slightly Xs/X —0.055. 

VI. SUMMARY 

Oblique plasmoid instabilities were analyzed within 
the context of the reduced MHD approximation. This 
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FIG. 8. (Color online) Three-dimensional flux ropes obtained from isosurfaces of the flux function, along with corresponding 
magnetic islands obtained from 2D cuts in the x — y plane. The x direction is vertical and z direction mostly into the page in 
the figure. Here, = 10^, Sl = 10^, and Bpo/B^o = 0.1. The top row is a parallel mode, where ^ = 0, and the bottom row is 
an oblique mode, where = 0.06 radians. The colorbar corresponds to colors in the 2D plots, and represents values of constant 
flux normalized by BpoX. 



required accounting for the 3D effect of wave variations 
in the guide field direction, which is the ignorable direc- 
tion in the 2D theory. An important difference between 
the 2D and 3D theories is the location of the resonant 
surfaces (where F = k • = 0). Considering a Harris 
equilibrium with guide field, this can take any value in 
3D: Xs = —XdiTctdinh[kzBzo/{kyBpo)]. In the 2D case, 
kz = and the resonant surface is always the null sur- 
face of the poloidal equilibrium field. The same is true 
in the 3D case only if there is no guide field present. The 
boundary layer analysis changed primarily in the outer 
region, where A' depends on the angle of obliquity [see 
Eq. (11)]. The approximate A' expression we obtained 
for a Harris equilibrium was shown to compare favorably 
with a numerical solution of the ideal MHD force balance. 
In the inner region, the angle of obliquity only entered in 
locating the resonant surface when evaluating F\xs). 

We found that unstable modes are confined to small 
angles \0\ ^ arctan(5po/52o). In the constant-?/^ regime, 
the most unstab le tearing m ode is an oblique mode 
\0\ {Bpo/Bzo)V{l^k^X^)/3 / 0. In the nonconstant- 

regime, the most unstable tearing mode is a parallel 
mode (^ = 0). The most unstable wavenumber is located 
at the intersection of these two regimes. The growth 
rate for this wavenumber is largest in the parallel direc- 
tion. By choosing an appropriate wavenumber at fixed 
Lundquist number (or vice versa) it was shown that a 
situation can arise in which only oblique modes are un- 



stable. The boundary layer theory was shown to compare 
well with numerical solutions of the linear reduced MHD 
equations. 
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